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Abstract 

The production of the neutral strange hadrons Kg, A and A has been measured in 
ep collisions at HERA using the ZEUS detector. Cross sections, baryon-to-meson 
ratios, relative yields of strange and charged light hadrons, A (A) asymmetry 
and polarization have been measured in three kinematic regions: Q 2 > 25 GeV 2 ; 
5 < Q 2 < 25 GeV 2 ; and in photoproduction (Q 2 ~ 0). In photoproduction the 
presence of two hadronic jets, each with at least 5 GeV transverse energy, was 
required. The measurements agree in general with Monte Carlo models and are 
consistent with measurements made at e + e~ colliders, except for an enhancement 
of baryon relative to meson production in photoproduction. 
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1 Introduction 



Production of K s , A and A hadrons has been extensively studied at particle colliders: 
e + e~ [1-3], ep [4], pp [5,6] and pp [7]. The data have been used to test QCD and build 
phenomenological models extending QCD predictions beyond what can be calculated from 
first principles. 

The results on Kg, A, and A production presented in this paper are based on a data 
sample of 121 pb -1 collected by the ZEUS experiment at HERA, about 100 times larger 
than used in previous HERA publications [4] and extend the kinematical region of the 
measurements, thereby providing a tighter constraint on models. 

The measurements have been performed in three different regions of Q 2 , where Q 2 is the 
virtuality of the exchanged boson: Deep Inelastic Scattering (DIS) with Q 2 > 25 GeV 2 ; 
DIS with 5 < Q 2 < 25 GeV 2 ; and photoproduction, Q 2 ~ GeV 2 , in which a quasi- real 
photon interacts with the proton. In the photoproduction sample, two jets, each of at 
least 5 GeV transverse energy, were required. 

The following measurements are presented in this paper: differential cross sections, baryon- 
antibaryon asymmetry, baryon-to-meson ratio, ratio of strange-to-light hadrons, and the 
A and A transverse spin polarization. There was no attempt to separate direct production 
from resonance decays: all sources for K s , A, and A production were included. These 
measurements are relevant for modeling production of hadrons at high energies, for ex- 
ample in Monte Carlo (MC) programs [8-12], and for testing the mechanism for baryon 
transport along the rapidity axis [13], the mechanisms for baryon production [14], effects 
due to QCD instantons [15-17] and the mechanisms for the transverse spin polarization 
of hadrons [18-20]. 



2 Experimental setup 

The data were collected by the ZEUS detector at the HERA ep collider during the running 
period 1996-2000. The data correspond to an integrated luminosity of 121 pb" 1 , of which 
82 pb" 1 were collected at y/s = 318 GeV (the electron or positrorij beam energy, E^ e£Lm , 
was 27.5 GeV and the proton beam energy was 920 GeV) and 39 pb -1 at s/s = 300 GeV 
(where the proton beam energy was 820 GeV). 

A detailed description of the ZEUS detector can be found elsewhere [21]. A brief outline 
of the components that are most relevant for this analysis is given below. 

1 In the text, electron beam, as well as scattered electron, apply to both electron and positron. 
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Charged particles are tracked in the central tracking detector (CTD) [22], which oper- 
ates in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD 
consists of 72 cylindrical drift chamber layers, organized in 9 superlayers covering the 
polar-angl^l region 15° < 6 < 164°. The transverse-momentum resolution for full-length 
tracks is a(p T )/p T = 0.0058p T © 0.0065 © 0.0014/pr, with p T in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [23] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections. The 
smallest subdivision of the calorimeter is called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, are a(E)/E = 0.18 /VE for electrons and 
a{E)/E = 0.35/ \[E for hadrons, with E in GeV. 

A three-level trigger system was used to select events on-line [21]. At the third level, DIS 
events were accepted on the basis of the identification of a scattered electron candidate 
using localised energy deposition in the CAL. As there was no possibility to select inclusive 
photoproduction sample, the requirement for photoproduction events [24] was based on 
running a jet algorithm using the energies and positions of the CAL cells. Events with 
at least two jets were accepted, where each jet was required to have transverse energy 
greater than 4.5 GeV and pseudorapidity less than 2.5. 

The luminosity was measured using the bremsstrahlung process ep — > ep^f with the lumi- 
nosity monitor [25], a lead-scintillator calorimeter placed in the HERA tunnel at Z = — 107 
m. 

3 Event reconstruction and selection 
3.1 Deep inelastic scattering sample 

The DIS events are characterised by a scattered electron detected in the CAL. The scat- 
tered electron was identified from the energy deposit in the CAL using a neural net- 
work [26]. The Bjorken variable xbj [27] and Q 2 were reconstructed using the double 
angle method (DA) [28] which has the best resolution in the Q 2 region studied. The 
inelasticity variable, y, was reconstructed using both the electron (e) [28] and Jacquet- 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the center of HERA. The coordinate origin is at the nominal interaction point. 

3 The pseudorapidity rj is defined as r\ = — lntan(0/2), where 9 is a scattering angle. 
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Blondel (JB) [29] methods. The following requirements were applied offline to select 
events with Q 2 > 25 GeV 2 (called the high-Q 2 sample): 

• | Z vtx | < 50 cm to reduce the background from non-ep collisions; 

• 38 < 5 < 65 GeV, where 5 = 5ZiC^» — Pz,i) and the sum runs over the energy and 
longitudinal momentum of all CAL cell deposits. This cut reduced the background 
from photoproduction and events with large radiative corrections; 

• an identified scattered electron with energy above 10 GeV; 

• the impact position of the scattered electron on the CAL satisfied \JX 2 + Y 2 > 36 cm; 

• the electron was isolated: the energy from all CAL cell deposits not associated with 
the scattered electron but in an 77 — cone of radius 0.8 centered on the electron was 
required to be below 5 GeV. This requirement reduced photoproduction background; 

• a track match with any electron falling in the range 0.3 < 9 < 2.6, well within the 
CTD acceptance. For 9 outside this region, 8 > 44 GeV was required. This cut further 
suppressed events from non-ep interaction and photoproduction; 

• Vjb > 0.02 to improve the accuracy of the DA reconstruction; 

• y e < 0.95 to remove events where fake electrons were found in the FCAL; 

• Ql A > 25 GeV 2 . 

The same selection was used to obtain the low-Q 2 DIS sample, except for the Q 2 require- 
ment and the position of the scattered electron, which were as follows: 

• the impact position of the scattered electron on the CAL was required to be outside 
a rectangle of dimensions 26 x 14 cm 2 , centred on the beam pipe; 

• 5 < Q 2 DA < 25 GeV 2 . 

The trigger for selecting low-Q 2 events was normally prescaled, so the data correspond to 
an integrated luminosity of 16.6 pb~ x . 

It should be noted that there was no jet requirement in the DIS event samples. 
3.2 Photoproduction sample 

Photoproduction events were selected applying the following criteria, described in an 
earlier publication [24]: 

• \Z vtx \ < 50 cm, to reduce background from non-ep collisions; 

• Vjb > 0.2, to further reduce background from non-ep collisions; 
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• events were removed where an electron was found with y c < 0.85, reducing the back- 
ground from neutral current DIS events; 

• Z/jb < 0.85, to reduce background from neutral current DIS events where the electron 
was not identified; 

• charged current DIS events were removed by rejecting events with P™ 1SS / y/Er > 
2.0 GeV 5 , where P™ 1SS is the missing transverse momentum and E? is the total trans- 
verse energy. 

Energy Flow Objects, reconstructed from a combination of calorimeter and tracking infor- 
mation to give the best resolution of kinematic variables, were used as the input to the kr 
cluster jet-finding algorithm [30], which was run in the longitudinally invariant inclusive 
mode [31]. The transverse energy of the jets was corrected for energy losses in inactive 
material in front of the CAL, as described in a previous publication [32]. An event was 
accepted if it contained at least two jets, both satisfying the following criteria: 

• jet transverse energy Ej, > 5 GeV; 

• jet pseudorapidity \r]i ct \ < 2.4. 

Photoproduction events selected in this way contributed about 10% to the total photo- 
production cross section. The photoproduction sample was divided into subsamples using 
the variable x® BS , a measure of the fraction of the photon energy transferred to the dijet 
system, defined as: 

x OBS _ Z^^T e 

7 22/ JB £ e beam 

where the sum runs over the two jets with highest transverse energy. In leading-order 
QCD, x° BS = 1 corresponds to direct photon processes in which the photon takes part 
in the hard scattering as a point-like particle. Resolved photon processes, in which the 
photon acts as a source of partons, populate the region at x^ BS < 1. The sample with 
x® BS > 0.75 is classified as direct-enriched, and that with x° BS < 0.75 as resolved- 
enriched. 



3.3 Strange particle reconstruction 

Candidates for long-lived neutral strange hadrons decaying to two charged particles were 
identified by selecting pairs of oppositely charged tracks, fitted to a displaced secondary 
vertex. Events were required to have at least one such candidate. These secondary 
vertices were found by the ZEUS track-finding software, which is based on minimizing 
the x 2 arising from fitting tracks to vertices [33]. Displaced vertices were typically more 
than 3 cm away from the primary vertex. The minimal distance required to resolve a 
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displaced vertex from the primary vertex was about 1 cm. The tracks fitted to this vertex 
were required to pass through at least the fifth superlayer of the CTD, the transverse 
momentum was required to be greater than 150 MeV and the absolute pseudorapidity 
in the laboratory frame was required to be less than 1.5. These constraints ensured a 
good track resolution and acceptance. The A, A and K° s particles may also be created in 
interactions with the beam pipe. To remove these events a collinearity cut on the angle 
between the reconstructed candidate momentum and the vector joining the primary vertex 
to secondary vertex was applied. This angle was restricted to be less than 0.2. 

The A(A) candidates were reconstructed by their charged decay mode to pvr"( pn + ) 
(branching ratio 63.9 ± 0.5% [34]). The track with the larger momentum was assigned the 
mass of the proton, while the other was assigned the mass of the charged pion, as the decay 
proton always has a larger momentum than the pion, provided the A(A) momentum is 
greater than 0.3 GeV. Additional requirements to select A(A) are given in the following: 

• 0.6 < P^ AB (A, A) < 2.5 GeV, where P^ AB (A, A) is the transverse momentum of the 
reconstructed candidate; 

• | rj LAB (A, A) |< 1.2, where ?7 LAB (A, A) is the pseudorapidity of the reconstructed 
candidate in the laboratory; 

• M(e + e~) > 0.05 GeV, to eliminate electron pairs from photon conversions 0; 

• M(iT + 7r~) < 0.475 GeV, to remove K$ contamination; 

• 1.11 < M(pn) < 1.122 GeV. 

The Kg meson candidates were reconstructed from the decays to 7r + 7r~ (branching ratio 
68.95 ± 0.14% [34]). Both tracks were assigned the mass of the charged pion. Additional 
requirements to select Kg are given in the following: 

• 0.6 < P^ AB (Kg) < 2.5 GeV, where P^ AB (Kg) is the transverse momentum of the 
reconstructed candidate; 

• | r] LAB (Kg) |< 1.2, where r] LAB (Kg) is the pseudorapidity of the reconstructed candi- 
date in the laboratory; 

• M(e + e~) > 0.05 GeV; 

• M(p7i) > 1.125 GeV, to remove A and A contamination. Here the mass of the proton 
was assigned to the track with larger momentum and the mass of the pion to the other 
track; 

• 0.48 < M(tx + ix~) < 0.52 GeV. 

4 M(ab) is defined as the invariant mass for two vertex tracks with the assignment of masses of particles 
a and b. 
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The mass peaks for Kg and A + A in the high-Q 2 sample are shown in Fig. [TJ The decay 
of Kg, A, and A was well understood as can be demonstrated by Fig. [21 which shows the 
proper decay times, reconstructed from the three dimensional decay length, compared to 
the expectations from MC simulation (see below). 

4 Event simulation 

Production of K s , A and A hadrons was modelled using the MC programs described 
below. In these models, strange quarks can be produced perturbatively by the boson-gluon 
fusion process (jg — > ss) or by gluon splitting in so-called parton showers. They may 
also originate from the proton parton densities or can be generated in non-perturbative 
string fragmentation. Strange hadrons are produced during hadronization, when quarks 
recombine into hadrons, and through the decays of other hadrons. Samples of events were 
generated to determine the response of the detector and obtain the correction factors 
required to convert the detector-level distributions to the hadron level. The generated 
events were passed through a full simulation of the detector, using Geant 3.13 [35], and 
processed with the same reconstruction program as used for the data. 

The high-Q 2 and the \ow-Q 2 DIS data were corrected to the hadron level using the Ari- 
adne 4 [8] MC program interfaced to Heracles 4.6.1 [36-39] via Djangoh 1.1 [39], to 
include QED corrections. The CTEQ proton parton density functions were used [40,41]. 
Ariadne is based on the Colour Dipole Model in which most QCD coherence effects are 
modelled as gluon emission from colour dipoles between partons. The program uses the 
Lund string model [42] to simulate the fragmentation of the partons. A significant param- 
eter governing the production of strange hadrons is the strangeness-suppression factor, X s , 
that is probability to produce s-quark pairs relative to u- and <i-quark pairs in the string 
fragmentation. This was set to 0.3, the default value, as found in e + e~ annihilation [43]. 
Other parameters that control baryon production were set to their default values^ [44]. 
This Ariadne sample was also used to compare to the final cross sections and ratios. 

A value of A s smaller than 0.3 is often preferred [2,4,45] for K s production. Therefore, 
a further sample with X s = 0.22 [46] was also generated and used for comparison. The 
DIS data were additionally simulated using the Lepto 6.5 MC program [9], which is 
based on first-order matrix elements plus parton showers (MEPS). The same Lund string 

5 The key parameters for baryon production in JETSET are the diquark-antidiquark pair production 
suppression PARJ(1)=0.10, the suppression of s quark pair production compared to u or d pair pro- 
duction X s = PARJ(2)=0.30, the extra suppression of strange diquark production compared with the 
normal suppression of strange quarks PARJ(3)=0.4 and the suppression of spin 1 diquarks compared 
with spin ones PARJ(4)=0.05. 
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model was used for the hadronisation, with X s = 0.3, and the same proton parton density 
functions as in the Ariadne sample. This was used for further comparison to the data. 

The photoproduction data were corrected to the hadron level using the Pythia 6 event 
generator [11], which consists of leading-order matrix-element calculations with initial- 
and final-state parton showering to simulate higher-order processes. The proton and 
photon PDFs were taken from GRV [47] and SaS2D [48] respectively. Multiple interac- 
tions [24,49,50], where more than one pair of partons (one parton from the photon and one 
parton from the proton) interact independently, were included. The default implementa- 
tion was used, with the p™ m [44] value of 2.7 GeV. The hadronisation is performed by the 
Lund string model, as in Ariadne, with the same parameters controlling the production 
of strange hadrons. Direct and resolved events were generated separately. For correction 
of the data, the direct and resolved subsamples were combined such that they gave a best 
fit to the data x° BS distribution. The Pythia sample was also used to compare to the 
final cross sections and ratios, in which case the direct and resolved events were combined 
according to the predicted cross sections. 



5 Cross-section determination 

The cross sections in the high-Q 2 DIS sample were measured in the kinematic region Q 2 > 
25 GeV 2 and 0.02 < y < 0.95. The cross sections in the low-Q 2 DIS sample were measured 
in the kinematic region 5 < Q 2 < 25 GeV 2 and 0.02 < y < 0.95. The cross sections in 
the photoproduction sample were measured in the kinematic region Q 2 < 1 GeV 2 and 
0.2 < y < 0.85, with the additional requirement of 2 jets, both satisfying Ej, > 5 GeV 
and |?7j et | < 2.4. In all three samples there was a further kinematic requirement that 
0.6 < P^ AB (K° S , A, A) < 2.5 GeV and \r] LAB (K s , A, A)| < 1.2. In all samples the measured 
cross sections were the luminosity-weighted average of the cross sections at the centre- 
of-mass energies yfs = 318 GeV and \fs = 300 GeV. In the DIS samples differential 
cross sections were measured as functions of P^ AB (K S , A, A), r] hAB (K s , A, A), ccbj, and 
Q 2 . In the photoproduction sample differential cross sections were measured as functions 
of P^ AB {K° S , A, A), V LAB {K S , A, A), and x° BS . 

The Kg, A and A differential cross sections in any variable Y were calculated using a 
standard bin-by-bin correction as follows: 

da _ N 

dY ~ A-C-B-AY' 

where N is the number of Kg, A or A in a bin of width AY, £ is the luminosity, A is the 
acceptance and B is the branching ratio. The acceptance was calculated from the MC 
samples described in Section HI and took into account migration effects and efficiencies 
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for each bin. The acceptance for each particle was calculated for each bin of a 6 x 6 grid 
in px and 77. An additional acceptance correction was made to obtain the differential 
cross sections as functions of variables other than p? or 77. In each bin a mass-sideband 
subtraction method was used to subtract the remaining combinatorial background in the 
Kg, A and A samples, which was at the level of ~ 3% in the Kg sample and ~ 6% in the 
A and A samples. 

The differential cross sections, j^, for the considered particle as a function of 77 and Pt 
can be converted into the multiplicity distribution j^, where a inc is the inclusive 
event cross section using the following factors: 0^ = 2.6 x 10~ 5 pb~ 1 for the high-Q 2 
sample, 8.9 x 10 -6 pb _1 for the low-Q 2 sample and 1.6 x 10~ 5 pb _1 for the photoproduction 
sample. On average, in the measured Pt and 77 region, for the high-Q 2 sample, there were 
about 0.017 A (or A) and 0.09 Kg per event. The corresponding multiplicities for the 
photoproduction sample were 0.077 for A (or A) and 0.27 for Kg. 

5.1 Systematic uncertainties 

The main sources which contributed to the systematic uncertainty were investigated by 
changing the analysis procedure, as outlined below, and observing the difference from the 
primary result. The total systematic uncertainty in the DIS samples for each bin was cal- 
culated by adding the individual contributions from the different variations in quadrature. 
The following sources of systematic uncertainties were considered: the energy measure- 
ment of the scattered electron; the measurement of Q 2 , 5, Z vtx and y; measurements of 
secondary-track pt and 77; the removal of A and Kg due to the collinearity cut; the impact 
position of the electron on the CAL; the estimation of the background. The systematic 
uncertainty on the cross sections was generally < 5%. Additionally, a 2% uncertainty due 
to the luminosity measurement was included for the cross sections. The branching-ratio 
uncertainties were deemed negligible and not taken into account. 

The most significant systematic error in the photoproduction sample was due to the 
uncertainty in the calorimeter energy scale (±3%), shown as a separate band in the 
figures. Typical uncertainties in the cross sections were < 10%, except in the highest bin 
of x® BS where the uncertainty was up to 80%. A systematic uncertainty of 7% on the 
photoproduction cross sections due to trigger efficiencies was included. The limitations of 
the Pythia Monte Carlo simulation, particularly in describing the x® BS distribution in 
the data, introduced a possible additional systematic uncertainty in correcting the data 
from detector level to hadron level. The effect of this was investigated by reweighting the 
Pythia x° bs distribution to the data and using this reweighted sample to correct the 
data. The difference between the results obtained with the reweighted sample compared 
to the primary results was treated as a further systematic uncertainty; uncertainties in 
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the cross sections were < 5%. 

The uncertainty of the tracking simulation [51] was negligible compared to all other 
sources. Most of the uncertainties discussed above cancel in the ratios and asymmetries 
presented in this paper. 

6 Results and discussion 

6.1 Cross sections 

Measured differential cross sections for the production of K s and A + A are shown in 
Figs. [3] - E] for the DIS data and in Figs. [7] and M for the photoproduction sample. The 
DIS cross sections are compared to the absolute predictions of Ariadne and Lepto MC 
calculations. The photoproduction cross sections are compared to the prediction of the 
Pythia MC with multiple interactions normalised to the data cross section. 

The Ariadne program with strangeness-suppression factor of 0.3 describes the data rea- 
sonably well, although the K s cross section for the high-Q 2 sample is overestimated. The 
slope of the P^ AB dependence is incorrect and the cross section at low XBj is underesti- 
mated for both low- and high-Q 2 samples. A similar comment can be made for the A + 
A cross sections. 

The description of the data by Ariadne with the strangeness-suppression factor of 0.22 
is less satisfactory but, as X s is not the only parameter influencing the cross section, a 
conclusion on what value of the X s can describe the data best can not be drawn. The 
Lepto MC does not describe the data well. 

In photoproduction, Pythia with multiple interactions describes the shape of the data 
dependence on P^ AB and r/ LAB adequately, but the x® BS dependence is predicted to be 
too flat and at the smallest x^ BS , in the resolved photon region, the description is poor. 

6.2 Baryon-antibaryon asymmetry 

The baryon-antibaryon asymmetry A is defined as: 

_ N(A) - N(A) 
N^ + NiA)' 

where A r (A), N(A) is the number of A and A baryons, respectively. 

The baryon-antibaryon asymmetry A has been measured and compared to MC predic- 
tions from Ariadne and Pythia (with X s = 0.3 for both, and also with A s = 0.22 for 
Ariadne). The following values were obtained: 
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• at high Q 2 : ^4 = 0.3 ± 1.3+°;g%, compared to the Ariadne (X s = 0.3) prediction of 
0.4 ±0.2%; 

• at low Q 2 : A = 1.2 ± compared to the Ariadne (A s = 0.3) prediction of 
1.0 ±0.2%; 

• in photoproduction: A = —0.07 ± 0.6^}'q%, compared to the Pythia prediction of 
0.6 ±0.1%. 

Figures [9] and [10] show the baryon-antibaryon asymmetry at high-Q 2 as a function of 
P^ AB , ?7 LAB , XBj, Q 2 and in photoproduction as a function of P^ AB , r/ LAB and x° BS . In all 
cases, the average baryon-antibaryon asymmetry is consistent both with no asymmetry 
and consistent with the very small asymmetry predicted by Monte Carlo. This suggests 
that in the considered parts of the ep phase-space, to a good approximation, baryons and 
antibaryons are produced according to the same mechanism. 

A positive asymmetry of 3.5% is predicted in DIS [13], due to the so called gluon-j unction 
mechanism that makes it possible for the "baryon number to travel" several units of 
rapidity, in this case from the proton beam direction to the rapidity around in the 
laboratory frame in which the measurements were made. Such a class of models can 
describe the significant positive baryon-antibaryon asymmetry which has been measured 
at the heavy-ion collider RHIC [52]. Although this prediction is not Q 2 dependent, it 
is not clear whether it could be extended down to Q 2 = and applied to the selected 
photoproduction sample. Combining statistical and systematic errors from the three 
samples, the average asymmetry is in disagreement with the 3.5% value. 

Although Ariadne predicts that on average about 15% (40% at the highest XBj for the 
low-Q 2 sample) of events with a reconstructed A or A originate from the exchanged photon 
coupling to an s or s quark from the proton, this measurement, being at low XBj, is not 
sensitive to the strange-quark asymmetry in the proton structure function as studied by 
NuTeV [53]. 

As our baryon-antibaryon asymmetry is consistent with no asymmetry, the A and A 
samples were combined together (except for the transverse polarization measurement) 
and results presented for the combined sample. 

6.3 Baryon-to-meson ratio 

The baryon-to-meson ratio TZ is defined as: 

_ A(A) + A(A) 
N(K S ) ' 

where N(A), N(A), N(K S ) refer to the number of indicated hadrons. 
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Figures [TT1 and [T2l show the measured and predicted TZ for the DIS samples as a function of 
P^ AB , r] LAB , xbj and Q 2 . The inaccuracies in describing the Kg and A + A cross sections 
for high Q 2 , mentioned earlier, are clearly reflected in TZ, but overall the Ariadne MC 
with X s = 0.3 follows the shape of the data distributions and is usually in agreement 
to better than about 10%. The low-Q 2 sample is described by the same Monte Carlo 
programs with even better accuracy. In order to have a better understanding of how 
TZ depends on xbj and Q 2 , TZ is shown as a function of Q 2 for fixed bins in xbj and as 
a function of xb] for fixed bins in Q 2 in Figs. [13] and [14j A dependence on Q 2 and a 
discrepancy between the data and MC can now be seen for the bins of higher xbj- For 
the two bins of higher Q 2 , the MC underestimates the data at low xej by up to 20%. 
The 1Z value varies between about 0.2 and 0.5, and is about 0.4 to 0.5 at low and 
low Q 2 . These values can be compared to measurements at e + e~ colliders, where for 
centre-of-mass energies from 10 to 200 GeV, TZ varies between about 0.2 and 0.4 [34]. 

Figure [151 shows TZ for the photoproduction sample. For the direct-enriched sample, where 
x° BS > 0.75, TZ is about 0.4, the same value as in DIS at low xbj and low Q 2 . However, TZ 
rises to a value of about 0.7 towards low x° BS (resolved-enriched sample), while it stays 
flat in the Pythia prediction. 

In order to study this effect further, the photoproduction events were divided into two 
samples. In the first, the jet with the highest transverse energy was required to contribute 
at most 30% to the total hadronic transverse energy. In this sample the events have largely 
isotropic transverse energy flow and therefore the sample is called "fireball-enriched" . The 
other sample, containing all the other events, was called "fireball-depleted" . Figure [ToTa) 
shows the distribution of events as a function of the total transverse energy, which is on 
average 30 GeV, and of the transverse energy of the jet with the highest . The line 
represents the cut used to separate fireball-enriched and fireball-depleted samples. Figure 
[ToTb) illustrates the fireball selection in relation to the fraction of the total transverse 
energy carried by two jets of the highest transverse energy. It can be seen that fireball- 
depleted events are dominated by dijet events carrying most of the total transverse energy 
and that the fireball-enriched and fireball-depleted samples have about the same number 
of events. 

The baryon-to-meson ratios for the fireball-enriched and fireball-depleted samples are 
presented in Fig. [17] for the data and Pythia MC. The measured TZ is larger for the 
fireball-enriched sample, most significantly at high P^ AB , than it is for the fireball-depleted 
sample. This feature is not reproduced by Pythia, which predicts almost the same TZ 
for both samples. The Pythia prediction reasonably describes the measured values of TZ 
for the fireball-depleted sample. This is not surprising as Pythia generates jets in events 
according to the multiple interaction mechanism [49], which makes several independent 
jets, like those in DIS or e + e~ where baryons and mesons are created locally. Provided 



11 



there is enough energy available, 1Z will be the same, regardless of the number of jets 
(ignoring some differences in quark and gluon fragmentation). 

Large values of TZ, larger than 1, have been measured at hadron and heavy- ion colliders: 
pp [6], pp [7,54] and RHIC [55-58]. 

6.4 Ratio of strange-to-light hadrons 

The ratio of strange-to-light hadrons was measured in terms of T: 

r _ N(K°s) 

where N(K S ) is the number of Kg and iV c h is the number of charged pions, charged kaons, 
protons and antiprotons, (excluding products of K s , A, and A decays) in the same region 
of P£ AB and r] LAB as the K° s . 

In Figs. dH] and [TH1 T is shown as a function of P^ AB and r] LAB for the high-Q 2 sample 
(for the low-Q 2 sample, not shown, the values are similar) and for the photoproduction 
sample. The MC predictions from Ariadne and Pythia are also shown. They follow 
the data reasonably well, preferring the strangeness-suppression factor to be smaller than 
0.3. The measured T lies between 0.05 and 0.1, varying with P^ AB for both the DIS and 
photoproduction. Similar values have been measured at e + e~ [34] for the ratio of the 
number of K s to the number of charged pions and are on average about 0.07 at centre- 
of-mass energies from 10 to 35 GeV, about 0.06 at the Z° and about 0.05 at 200 GeV. It 
can be concluded that T is about the same in e + e~ and ep. 

In order to see whether T depends on the transverse energy flow, the fireball selection, as 
discussed above, was applied to the photoproduction events. Figures flWc) and (d) show 
the measured and predicted T for the fireball selection. The quantity T hardly depends 
on the fireball selection, as predicted by Pythia. 

Fireball events are candidates for events where QCD instantons [15-17] could play a 
role [59-61], since they are characterised by isotropic transverse energy flow. Another 
expectation is a likely enhancement of heavier-quark production relative to light quarks, 
due to the required flavour democracy. Searches for QCD instantons in DIS have been 
reported by HI [62] and ZEUS [63]. No effect was identified due to QCD instantons, as 
the expected effects were small compared to the background at the relatively large Q 2 
required. Bigger effects are expected at lower Q 2 . 

If QCD instantons contribute to the fireball event sample, then T would be expected to 
be different, possibly larger, for the fireball events than the typical value of about 0.07 or 
0.08. As this is not the case, this measurement of T in photoproduction does not support 
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the idea that QCD instantons contribute significantly to the production of the fireball 
events. It should be noted that there is only a qualitative prediction on the contribution 
from QCD instantons based on democracy of flavours, including heavy flavours, subject 
to available energy. The only existing calculation [61] applies to DIS and only considers 
three massless flavours. Since there is no charm-quark contribution, this calculation is 
probably only valid at low particle multiplicities, where the number of Kg is predicted to 
be about twice as big as that predicted by Ariadne. 



6.5 Polarization 



In analogy with QED, the spin-orbit interaction leads to polarization in scattering due to 
the strong interaction [18-20]. Unpolarized s quarks get partially transversely polarized 
due to elastic scattering in the coloured field along the direction of hi x kf, where hi and 
kf stand for the initial and final momenta of the s quarks. The degree of the polarization 
depends on the scattering angle and the strength of the coloured field. In the constituent 
quark model, the A inherits its spin from the s quark, and kf is along the A momentum. 
As hi is unknown in our case, the electron beam direction was used instead (the effect of 
using the jet direction was also investigated). 

The transverse polarization V A (V A ) is defined by the form of the proton (antiproton) 
angular distribution: 

1 dN 1 



Ndcos# 2 
1 dN 1 



-[1 + aV A cos 6} 



[1 - aP A cosfl] 



Ndcosfl 2 



where a is the decay asymmetry parameter, measured to be a = 0.642 ± 0.013 [34], and 
9 is the angle between the proton (antiproton) momentum boosted to the rest frame of 
the A (A) and the polarization axis, /c^ eam x k^. An example of the angular distribution 
of the proton's (antiproton's) momenta with respect to the polarization axis, boosted to 
the A(A) rest frame, is shown in Fig. [2DJ 

Fitted values of the transverse polarization V A and V A are presented in Table [U for high- 
and low-Q 2 DIS and for photoproduction. All values are compatible with no polarization. 
Also presented are the polarization values obtained by investigating the angular distri- 
bution of the higher-momentum ir from Kg decays, as a further test of any systematic 
detector effect. 
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7 Conclusions 



Measurements of Kg, A and A production have been made at HERA, using 121 pb" 1 of 
data collected with the ZEUS detector. The following conclusions have been obtained: 

1. in high- and \ow-Q 2 DIS, Ariadne reproduces gross features of the cross sections but 
shows discrepancies in detail. Overall, the strangeness suppression factor A s = 0.3 is 
preferred to A s = 0.22. Pythia, normalised to the data, describes the dependence 
of the photoproduction cross sections on P^ AB and rj LAB satisfactorily but fails to 
reproduce the x° BS dependence; 

2. the numbers of A and A produced are consistent with being equal; 

3. except for the resolved photon interactions, the measured ratio of baryons to mesons 
71, defined as: 

_ N(A) + N(A) 
N(K°) ' 

is in the range between 0.2 and 0.5, similar to measurements at e + e~ colliders. Ari- 
adne and Pythia follow the shapes of the data on the selected observables but in 
many cases fail quantitatively at the 10 to 20% level; 

4. in the resolved photon region, the ratio of baryons to mesons is large, significantly 
larger than measured in e + e~ interactions and significantly larger than the Pythia 
prediction; 

5. the ratio of strange-to-light hadrons measured in terms of T: 

N(K°s) 



T 



N, 



ch 



is compatible with measurements at other colliders and is described by Ariadne and 
Pythia for all investigated samples of events. For the DIS sample, the strangeness 
suppression factor A s = 0.22 is preferred to the default value of X s = 0.3. For the 
photoproduction sample, Pythia with X s = 0.3 overestimates the data, but describes 
the shape of the distributions. There is no indication of any unusual yield of strange- 
hadrons in the fireball-enriched sample, as would be qualitatively expected had there 
been a significant contribution from QCD instantons; 

6. No evidence has been found for non-zero transverse polarization in inclusive A or A 
production. 
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Polarization (%) 




High-g 2 DIS 


Low-g 2 DIS 


Photoproduction 


A 


-1.3±4:.3{stat.)to 8 {syst.) 


-4.0 ± 5.3(stat.)tii{syst.) 


-2A±2.2{stat.) 


A 


-2.2 ±4:.2(stat.)tl A 3 (syst.) 


-8.5±5.5(stat.)t 4 2 7 1 (syst.) 


-5.8 ± 2.2(stat.) 




-1.5 ± l.l(stat.) 


-0.05 ± 1.5{stat.) 


-0.5 ± 0.2(stat.) 



Table 1: The transverse polarization values for A and A, expressed here in 
%, in the high-Q 2 DIS ( Q 2 > 25 GeV 2 and 0.02 < y < 0.95/, low-Q 2 DIS ( 
5 < Q 2 < 25 GeV 2 and 0.02 < y < 0.95), and photoproduction ( Q 2 < 1 GeV 2 , 
0.2 < y < 0.85 and with two jets Ej, > 5 GeV and \rf ct \ < 2A) samples. Only A 
and A in the range 0.6 < P^ AB < 2.5 GeV and |r; LAB | < 1.2 are considered. The 
statistical error is quoted for all samples, together with the systematic uncertainty 
associated with the measurement for the high-Q 2 and low-Q 2 samples. A similar 
systematic uncertainty is expected for the photoproduction sample. Also shown, as 
a test of any systematic effect, are the polarization values obtained by investigating 
the angular distribution of the higher-momentum re from K® decays. 
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Figure 1: Mass distributions of the secondary vertex candidates in the A + A 
and Kg samples assuming pn and txtx decays respectively. Only candidates in the 
range 0.6 < P^ AB < 2.5 GeV and \t] lab \ < 1.2 for events with Q 2 > 25 GeV 2 and 
0.02 < y < 0.95 are displayed. Statistical errors are smaller than the point size. 
The number of total A + A and K s candidates located within the vertical lines are 
estimated to be 10731 and 73140, respectively, after the background subtraction. 
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Figure 2: Detector-level distributions of ct, where t is the reconstructed proper 
lifetime, for A, A and K s samples for data and Ariadne. The Ariadne histogram 
is normalised to the same number of events as the data. Statistical errors are shown 
unless smaller than the point size. 
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Figure 3: Differential Kg cross-sections as a function of and r] in the 

range 0.6 < P? AB < 2.5 GeV and \t] lab \ < 1.2 for events with 5 < Q 2 < 25 GeV 2 , 
0.02 < y < 0.95 and Q 2 > 25 GeV 2 , 0.02 < y < 0.95. Statistical errors (inner 
error bars) and the systematic uncertainties added in quadrature are shown, unless 
smaller than the point size. The histograms show predictions from Ariadne and 
Lepto using the stated strangeness suppression. 
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Figure 4: Differential Kg cross sections as a function o/xej and Q 2 in the range 
0.6 < P^ AB < 2.5 GeV and |r/ LAB | < 1.2 for events a) with 5 < Q 2 < 25 GeV 2 , 
0.02 < y < 0.95 5; Q 2 > 25 GeF 2 , 0.02 < y < 0.95 and c) Q 2 > 5 GeF 2 , 
0.02 < y < 0.95. Statistical errors and the systematic uncertainties added in 
quadrature are smaller than the point size. The histograms show predictions from 
Ariadne and Lepto using the stated strangeness suppression. 
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Figure 5: Differential A + A cross sections as a function of P^ AB andr] hAB in the 
range 0.6 < P? AB < 2.5 GeF and |?7 LAB | < 1.2 /or events with 5 < Q 2 < 25 GeF 2 , 
0.02 < y < 0.95 and Q 2 > 25 GeV 2 , 0.02 < y < 0.95. Statistical errors (inner 
error bars) and the systematic uncertainties added in quadrature are shown, unless 
smaller than the point size. The histograms show predictions from Ariadne and 
Lepto using the stated strangeness suppression. 
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Figure 6: Differential A+ A cross sections as a function of x-q] andQ 2 in the range 
0.6 < P^ AB < 2.5 GeV and |r/ LAB | < 1.2 for events a) with 5 < Q 2 < 25 GeV 2 , 
0.02 < y < 0.95 b) Q 2 > 25 GeV 2 , 0.02 < y < 0.95 and c) Q 2 > 5 GeF 2 ; 
0.02 < y < 0.95. Statistical errors and the systematic uncertainties added in 
quadrature are smaller than the point size. The histograms show predictions from 
Ariadne and Lepto using the stated strangeness suppression. 
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Figure 7: Differential Kg cross sections as a function of Pf , r t i and , 
in the range 0.6 < P^ AB < 2.5 GeV and \r] LAB \ < 1.2 for events with Q 2 < 1 GeV 2 , 
0.2 < y < 0.85 and at least two jets both satisfying > 5 GeV and \rf ct \ < 2.4. 
Statistical errors are shown, unless smaller than the point size, together with the 
systematic uncertainty arising from the trigger efficiency added in quadrature. The 
uncertainty arising from the jet energy scale is also shown (shaded band). The 
solid histogram shows the prediction from Pythia (with multiple interactions), 
normalised to the data. 
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Figure 8: Differential A + A cross sections as a function ofPji , r/ 



and x 



OBS 



in the range 0.6 < P? AB < 2.5 GeV and \rj LAB \ < 1.2 for events with Q 2 < 1 GeV 2 , 
0.2 < y < 0.85 and at least two jets both satisfying > 5 GeV and \r]^ ct \ < 2.4. 
Statistical errors are shown, unless smaller than the point size, together with the 
systematic uncertainty arising from the trigger efficiency added in quadrature. The 
uncertainty arising from the jet energy scale is also shown (shaded band). The 
solid histogram shows the prediction from Pythia (with multiple interactions), 
normalised to the data. 



28 



ZEUS 



~ i — I — i — i — i — i — I — i — i — i — r 
• ZEUS 96 -00 

LEPTO 0.3 

ARIADNE 0.3 

ARIADNE 0.2^ 




2 

_LAB 

P T (GeV) 



2.5 




0.25 

■ < 

Z 0.2 
+ 

z < 0.15 

T< 0.05 

\ 
5. -0.05 

-0.1 
-0.15 

-0.2 
-0.25 





0.25 






z 


0.2 


+ 




< 

z 


0.15 




0.1 




0.05 











'< 






-0.05 




-0.1 




-0.15 




-0.2 




-0.25 



Q 2 > 25 GeV 2 
0.02 < y < 0.95 



-1 



-0.5 



0.5 



1 

LAB 



10 



V Bj 



Q 2 (GeV 2 ) 



Figure 9: The ratio a function of P^ AB , r] LAB , XBj and Q 2 , in 

the range 0.6 < P^ AB < 2.5 GeV and |?7 LAB | < 1.2 for events with Q 2 > 25 GeV 2 
and 0.02 < y < 0.95. Statistical errors (inner error bars) and the systematic 
uncertainties added in quadrature are shown. The histograms show predictions 
from Ariadne and Lepto using the stated strangeness suppression. 
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0.2 < y < 0.85 and at least two jets both satisfying E^ 1 > 5 GeV and \if et \ < 
2.4. Statistical errors are shown, together with the uncertainty arising from the jet 
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Figure 11: The ratio n \k°) as a f unc ^ on °f Pt an< ^ rj in the range 
0.6 < P^ AB < 2.5 GeV and |r/ LAB | < 1.2 for events with 5 < Q 2 < 25 GeV 2 , 
0.02 < y < 0.95 and Q 2 > 25 GeV 2 , 0.02 < y < 0.95. Statistical errors (inner 
error bars) and the systematic uncertainties added in quadrature are shown. The 
histograms show predictions from Ariadne and Lepto using the stated strangeness 
suppression. 



31 



ZEUS 



z 

X 0.5 

■ < 

z 

+ . 0.4 

z 

0.3 





0.2 




0.1 









0.6 


z 










0.5 


z 




+ 

< 


0.4 


z 










0.3 




0.2 




0.1 








I I I I I lAl I 14. 




:(a) 



5 < Q 2 < 25 GeV 2 



0.5 

■ < 



0.4 
0.3 
0.2 
0.1 




ll I I I I I I I I I I I II I I 1 1 1 



(b) 



Q 2 > 25 GeV 2 



10 



10 



10 



10 



10 



V Bj 



V Bj 



Z ^1' 1^' Li 1 '4." "&"h/uJtwik 1 * _ - . j x 

' no " Ti " " " " " " *i n 1 1 1 1 1 - , £ , — x. 



4 I I I I II I I I H | | | | 



t t 

"I TTiiiiii, 



:(C) 



• ZEUS 121 pb" 1 
LEPTO (0.3) 

ARIADNE (0.3) 

ARIADNE (0.22) 



10 



10 



Q 2 (GeV 2 ) 



Figure 12: The ratio N ^^o^ as a function o/xbj and Q 2 in the range 0.6 < 
P^AB < 2 5 Ge y and ^LABi < l 2 j or evmts mfh a j 5 < q2 < 25 Q Q2 < 
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0.95. Statistical errors (inner error bars) and the systematic uncertainties added in 
quadrature are shown. The histograms show predictions from Ariadne and Lepto 
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Figure 13: The ratio N ^^o^ a function of Q 2 for four bins of XBj, in 
the range 0.6 < P^ AB < 2.5 GeV and \r] LAB \ < 1.2 for events with Q 2 > 5 GeV 2 



and 0.02 < y < 0.95. Statistical errors only are shown. The histograms show 
predictions from Ariadne with a strangeness-suppression factor of 0. 3. 
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and 0.02 < y < 0.95. Statistical errors only are shown. The histograms show 
predictions from Ariadne with a strangeness-suppression factor of 0. 3. 
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Figure 15: The ratio V(x u ) as a f unc ^ on °f Pj ' > r] and x® , in the 
range 0.6 < P^ AB < 2.5 GeV and \r] LAB \ < 1.2 for events with Q 2 < 1 GeV 2 , 
0.2 < y < 0.85 and at least two jets both satisfying > 5 GeV and \r]^ ct \ < 2.4. 
The statistical errors are shown, unless smaller than the point size. The shaded 
band shows the uncertainty arising from the jet energy scale. The solid histogram 
shows the prediction from Pythia (with multiple interactions). 
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Figure 16: (a) Distribution of photoproduction events as a function of the total 
transverse energy and the transverse energy of the jet with the highest transverse 
energy. Here the event sample was selected without any strange hadron require- 
ments. The solid line represents the cut used to separate fireball- enriched and 
fireball- depleted subsamples. (b) The fraction of the total Et carried by the two 
jets of highest transverse energy for the same data sample as in (a). Fireball- 
enriched (Ej , /Et < 0.3) and fireball-depleted (E^ /Et > 0.3) samples are 
shown. Statistical errors are smaller than the point size. 
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Figure 17: The ratio N ^^o^ a function of P^ AB , <q LAB and x® BS , in the 

range 0.6 < P^ AB < 2.5 GeV and \rj LAB \ < 1.2 for events with Q 2 < 1 GeV 2 , 0.2 < 
y < 0.85 and at least two jets both satisfying > 5 GeV and \rf ct \ < 2.4. The 
ratios from the fireball- enriched sample (squares) and the fireball- depleted sample 
(triangles) are shown for the data. The prediction from Pythia for the fireball- 
enriched (solid line) and for the fireball- depleted (dashed line) samples are shown. 
Statistical errors are shown. The highest x° BS bin (0.9 < x® BS < 1.0) of the fireball 
enriched sample is omitted due to insufficient statistics. 
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Figure 18: The ratio as a function of P^ Aa and ^ LAB ; in the range 

0.6 < P^ AB < 2.5 GeV and |r/ LAB | < 1.2 for events with Q 2 > 25 GeV 2 and 0.02 < 
y < 0.95. Statistical errors are smaller than the point size. The histograms show 
predictions from Ariadne and Lepto using the stated strangeness suppression. 
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Figure 19: The ratio N< ^ < ^ as a function of P? AB andrj LAB , in the range 0.6 < 
Pf AB < 2.5 GeV and \r] hAB \ < 1.2 for events with Q 2 < 1 GeV 2 , 0.2 < y < 0.85 
and at least two jets both satisfying > 5 GeV and |^ et | < 2.4. The ratio is 
shown for all events in (a) and (b) and for the fireball- enriched sample and the 
fireball- depleted sample in (c) and (d). Prediction from Pythia (with multiple 
interactions) predictions with a strangeness suppression factor of 0.3 are shown as 
solid and dashed histograms. Statistical errors are smaller than the symbols. 
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Figure 20: Angular distributions of the highest-momentum decay particle from 
a A or A in the range 0.6 < P^ AB < 2.5 GeV and |?7 LAB | < 1.2 for events with 
Q 2 > 25 GeV 2 and 0.02 < y < 0.95, where 9 is the angle between the decay- 
particle momentum vector and the polarization axis, in the rest frame of the A or 
A. Statistical errors (inner error bars) and the systematic uncertainties added in 
quadrature are shown. The first-order polynomial fit (solid line) from which the 
polarization is obtained is also shown. 
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